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Abstract 
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and overlap fermions as well as Wilson and C'(a)-improved Wilson fermions. 
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I. INTRODUCTION 



Lattice QCD calculations are necessarily performed with a finite lattice spacing, at a finite 
spatial extent and presently, with unphysically large quark masses. Thus, in order to obtain 
continuum, infinite volume results at the physical quark masses, multiple extrapolations 
are required. At large pion masses, naive (usually polynomial) extrapolations of results 
obtained with different quark masses and lattice spacings are performed. However, for results 
computed with quark masses in the chiral regime (corresponding to pion masses on the order 
of 400 MeV or lower), a generalization of chiral perturbation theory, which includes lattice 
spacing effects, can provide the analytic form of the dependence of many observables with 
the lattice spacing and quark mass [ll, HI, 0, 0, [sl, 0]. These formulae typically involve low- 
energy constants (LECs) parameterizing the physics of the QCD scale (and shorter) whose 
values are unknown. These LECs and corresponding operators (as well as quantum loop 
contributions) associated with lattice spacing effects, represent unphysical lattice spacing 
artifacts which must be removed from the correlation functions to obtain results in the 
continuum. The remaining LECs, which survive the continuum limit, represent the short- 
distance QCD physics of interest. One can view the role of the lattice calculations as 
determining the value of these physical constants, which when combined with the quantum 
loops, completely describe the low energy dynamics of QCD. 

The energy of a system can be obtained from the long time behavior of Euclidean cor- 
relation functions. In practice, baryonic observables (specially those connected to excited 
or multi-baryon states), suffer from a signal-to-noise degradation at large times while 



at short times, contamination from the many excited states masks the information about 
the low lying states. This naturally leads to the use of anisotropic lattice actions [sl, [9|, in 
which the lattice spacing in time is taken to be at = fls/C- Anisotropic lattices have been 
used extensively in the study of heavy quarks and quarkonia [l^, ll[ , glueballs [H, |l3[ , and 
excited state baryon spectroscopy 3. It is also anticipated that anisotropic lattices may 



aid in the study of nucleon-nucleon [15|, Il6[ and hyperon-nucleon [17| interactions on the 
lattice. For heavy systems made of light quarks, such as nucleon-nucleon systems, the rapid 
degradation of the signal-to-noise ratio of the correlation functions is particularly problem- 
atic in identifying the ground state. Thus, when the excited state contamination has finally 
died away, there are very few time slices left before the noise dominates the correlation func- 
tions, leaving only a very narrow window in time to study these systems. With the use of 
anisotropic actions, one may gain more useful information allowing an earlier identification 
of the ground state plateau. 

The production of a large ensemble of dynamical anisotropic Wilson gauge configurations 
is currently underway [l9|. It is therefore the aim of this work to construct the low energy 
effective theory which encodes these new anisotropic effects, so that they might be studied 
and systematically removed from the information extracted from the correlation functions 



generated with these lattices. We begin by constructing the Symanzik action [20|, |21[ for 



both the Wilson [22| and 0{a) improved Wilson 23| actions in Sec. [TTl Then in Sec. IIIIl 



we construct the chiral Lagrangians [2J, |25|, |26j relevant for these anisotropic lattices for 



both mesons and baryons, focussing on the new effects arising from the anisotropy. We also 



^ Recently, a "restless pions" boundary condition has been proposed [tJ designed to improve this signal- 
to-noise problem. Combined with anisotropic lattices, this may allow for a significant increase in the 
available time-slices to study the ground-state of these heavy systems. 
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provide extrapolation formulae for these hadrons with their modified dispersion relations. 
In Sec. IIII A H we highlight an important feature of anisotropic actions; for a fixed spatial 
lattice spacing and bare fermion mass, if the isotropic action is in the QCD-phase, this does 



not guarantee the anisotropic action is outside the Aoki-phase [27 
SeclYl 



We then conclude in 



II. ANISOTROPIC WILSON LATTICE ACTION 



The starting point for our discussion is the anisotropic lattice action and its symmetries, 



from which we will construct the continuum effective Symanzik action [20|, l2l|] , which will 
then allow us to construct the low-energy EFT describing the hadronic interactions inclu ding 
the dominant lattice spacing artifacts [2]. For the isotropic Wilson (and 0{a) improved 23 
Wilson) action, this program has been carried out to 0{a^) for the mesons [2[ B S] and 
baryons [s], [6|. This work is a generalization of the previous work, extending the low-energy 
Wilson EFT to include the dominant lattice artifacts associated with the anisotropy. 

The (!?(a)-improved anisotropic lattice action, in terms of dimensionless fields, is given 



by |28|,l29| 



t.i 



so 



i<3 



(1) 

(2) 
(3) 



Here, is the bare anisotropy, Pij and Pu are space-space and space-time plaquettes of the 
gauge links U . The bare (dimensionless) quark mass is atmo, Wt{U) and Wi{U) are Wilson 
lattice derivatives and u is a parameter which must be tuned to correct the "speed of light." 
The fields FuiU) and Fij{U) are lattice equivalents of the gauge field-strength tensor in 
the space-time and space-space directions. is the unimproved Wilson action and 5*1 
is the anisotropic generalization of the Sheikholeslami-Wohlert term. The coefficients, q 
and Cr, appearing in S*! are needed for 0{a) improvement of the anisotropic Wilson lattice 
action 23|. At the classical level, they have been determined to be |29|] 



Ct 



(4) 



where ^ = as/ at is the renormalized anisotropy.^ The choice of using ^ as opposed to is 
conventional and the difference amounts to a slightly different value of u. 

This anisotropic lattice theory retains all the symmetries of the Wilson action, except 
for the hypercubic invariance, and thus respects parity, time-reversal, translational invari- 
ance, charge-conjugation and cubic invariance. In addition, for suitably tuned bare fermion 



These parameters can also be tadpole improved with no more effort than in the isotropic action 3C | 
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masses, atruQ, the theory has an approximate chiral symmetry, SU{Nf)L ® SU{Nf)fi, which 
spontaneously breaks to the vector subgroup.^ 



A. Anisotropic Symanzik Action 



We begin by constructing the Symanzik Lagrangian for the unzmjorowed anisotropic Wilson 
Sl lattice action. This will allow us to set our conventions and introduce a new basis of 
improvement terms which is advantageous to studying the new lattice artifacts which are 
remnants of the anisotropy. In terms of dimensionful fields, the anisotropic Symanzik action 
is given by 

^Sym ~ J ^ ^^Sym ' 

'-Sym ~ '-Sym ^ "'s'-^Sym ' "'s'-Sym ' l^J 

where we have conventionally chosen to use the spatial lattice spacing as our Symanzik 
expansion parameter. In terms of the dimensionful fermion fields 

1 _ 1 - 

q ~ —ip , q ~ —i) , (6) 

the anisotropic Lagrangian is given through 0{a) by 

^\ym = q[Ip+mq]q + asq Ct auFu + c,. ^ crijFij q . (7) 

We have assumed that the parameter u has been tuned in such a way as to make the breaking 
of 0(4) symmetry to vanish in the continuum limit. Otherwise, the quark kinetic term would 
separate into two terms, with an additional free parameter appearing in eq. (171). To clearly 
identify the new lattice spacing effects associated with the anisotropy, it is useful to work 
with the basis, 

'^iym = q\Ip+mq] q + qcTf^^F^^ q + a^ict - c^) qauFu q , (8) 

from which we recognize the first 0{a) term as the Sheikholeslami-Wohlert [1^ term which 
survives the isotropic limit, csw = Cr/2. The second 0{a) term c|y^ = (q — Cr), is an 
artifact of the anisotropy and the focus of this work. We can classify the effects from this 
operator and subsequent anisotropy operators at higher orders into two categories: those 
which contribute to physical quantities in a fashion similar to the lattice spacing artifacts 
already present, and those which introduce new hypercubic breaking effects. The first type 
of effect will be difficult to distinguish from the already present lattice spacing artifacts 
which survive the isotropic limit. ^ The second category of effects are unique to anisotropic 
actions, and therefore more readily identifiable from correlation functions. 



^ As mentioned in the Introduction, in Sec. lIII A 1[ we will discuss Aoki-phases on anisotropic lattices which 
are not forbidden even if the equivalent isotropic action is not in the Aoki regime. 

^ Assuming a given set of lattice calculations are close to the continuum limit, and that a range of 
anisotropics is employed, one can disentangle the lattice artifacts associated with the new anisotropic 
operator from those which survive the isotropic limit. 
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A useful manner to quantify these new anisotropic effects is to recognize that the 
anisotropy introduces a direction into the theory, which we can denote with the four-vector 



n«=(l,0). (9) 



This allows us to re-write the anisotropic Lagrangian in a manner amenable to spurion 
analysis 

= q[p+mq]q + a^q Csw (y,,uF^,y + c^w'^i'^l (yf^xFyx q ■ (10) 

We then promote both UgCsw and asCg^ruj^ul to spurion fields, transforming under chiral 
transformations in such a way as to conserve chiral symmetry, 

dsCsw — ^ L{asCsw)R^ , {asCswY — ^ RiasCswYL^ (11) 

ascl^ului — > L{as4^ului)R^ , {ascl^yuluiy — > R{asc\w^lul)^ L'^ (12) 

By constraining agCsw and agC^sw'^lL'^v their hermitian conjugates to be proportional the 
flavor identity, they both explicitly break the SU{Nf)l ® SU{Nf)r chiral symmetry down 
to the vector subgroup, just as the quark mass term. In addition, we promote asc|^yM^M| 
to transform under hypercubic transformations, so as to conserve hypercubic symmetry, 

^sC^SW^ft^u ^ O-sC^t^WpIt^ A^pAj^o- . (13) 

By constraining = (1,0), this spurion explicitly breaks the hypercubic symmetry of the 
action down to the cubic sub-group. 

It is worth noting that, in fact, at 0{a), the anisotropic Symanzik action retains an 
accidental 0(3) symmetry in the spatial directions. Close to the continuum, isotropic lattice 
actions retain an accidental Euclidean 0(4) (Lorentz) symmetry as the operators required to 
break this symmetry are of higher dimension and thus become irrelevant in the continuum 



limit |20l. |21|. This phenomena is observed in the isotropic limit where the 0(4) symmetry is 
broken by the operator a?q'-)^D^D^D^q. For unimproved anisotropic Wilson fermions, the 
breaking of the the hypercubic to cubic symmetry (which can be viewed as the breaking of 
the accidental 0(4) to the accidental 0(3) symmetry) occurs one order lower in the lattice 
spacing, at 0(a), and therefore this will likely be a larger lattice artifact than the 0(4) 
breaking of the isotropic action. We now perform a similar analysis for the 0{a?) Symanzik 
action. 



1. 0{a?') Symanzik Lagrangian 

In Ref. [il, the complete set of 0{a?) operators in the isotropic Symanzik action for Wilson 
fermions was enumerated, including the quark bi-linears and four-quark operators. From an 
EFT point of view, it is useful to classify these operators in three categories, those operators 
which do not break any of the continuum (approximate) symmetries, those which explic- 
itly break chiral symmetry and those which break Lorentz symmetry. Most of the O(a^) 
operators belong to the first category. Because of their nature, they are the most difficult 
to determine and ultimately lead to a polynomial dependence in the lattice spacing of all 
correlation functions computed on the lattice (which can be parameterized as a polynomial 
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dependence in a of all the coefficients of the chiral Lagrangian). The second set of operators, 
those which explicitly break chiral symmetry, can be usefully parameterized within an EFT 
framework as is commonly done with chiral Lagrangians extended to include lattice spacing 
artifacts [2I, H, 0, HI, 0]. The last set of operators which break Lorentz symmetry, can also be 
usefully studied in an EFT framework. In the meson Lagrangian, these effects are expected 
to be small as they do not appear until 0{p^a^) while in the heavy baryon Lagrangian, 
these effects appear at 0{a?) j^. To distinguish these Lorentz breaking terms from the 
general lattice spacing artifacts appearing at O(a^), one must study the dispersion relations 
of the hadrons, and not merely their ground states. This is also generally true of all the 
anisotropic lattice artifacts which we now address. 

In the construction of the anisotropic action, it is also beneficial to categorize the operators 
into several categories along the lines of those in the isotropic action mentioned above. We 
do not show all of the new operators, as their explicit form will not be needed, but instead 
provide a representative set of the new anisotropic operators which illustrate the new lattice- 
spacing artifacts. In the ffist category, we begin with operators which in the isotropic limit 
do not break any of the continuum symmetries. Using the notation of Ref. [3], and using a 
superscript-'^ to denote the new operators due to the anisotropy, we find for example 

of ^ {Of , «0f } = {qD,IpD,q , qD.IpD.q} 

of ^ {Of , «of } = m,q){qi,q) , {qitq){qitq)} ■ (i4) 

In the second category, operators which explicitly break chiral symmetry, we find 

of {Of , «0f } = {{q<y,,q){qa,,q) , (g , (15) 

from which we observe that there is an operator which both breaks chiral and hypercubic 
symmetry. The last category of operators stems from the Lorentz breaking operator in the 
isotropic limit, 

of {Of , ^Of} = {q-f,D,D,D^ q , qi,D,D,D, q} , (16) 

from which we note that there is only one operator which breaks the accidental 0(3) sym- 
metry down to the cubic group, ^Of , and therefore the dominant 0(3) breaking artifacts in 
principle can be completely removed from the theory by studying the dispersion relation of 
only one hadron, for example the pion. Each of these new operators can be written in their 
spuriously hypercubic-invariant form by making use of the anisotropic vector we introduced 
in Eq. ©, 

^Of =ululqD,IpD,q, 
^Ofi = u^ui {qif,q){qiuq), 
^Of = u^ui {qaf,xq){q(ruxq), 

^Of = SlJlJlr6l,qi,'D^,D^,D,, g, (17) 

and similarly for the rest of the dimension-6 anisotropic operators, ^Of g, and ^Of_]^g. In 
this equation, for we have defined 

= V - ^I'^i ■ (18) 
Most of these operators do not break chiral symmetry, and therefore a re p resent for chirally 
symmetric fermions such as domain- wall 31, HI, 33 1 and overlap [13, 35, 3^ fermions. We 



now proceed to construct the anisotropic chiral Lagrangian. 
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III. ANISOTROPIC CHIRAL LAGRANGIAN 



Now that we have the complete set of Symanzik operators through O(a^) relevant to 
the anisotropic Wilson action and the 0{a) improved version thereof, we can construct the 
equivalent operators in the chiral Lagrangian which encode these new anisotropic artifacts. 
We begin with the meson Lagrangian and then move to the heavy baryon Lagrangian. 

A. Meson Chiral Lagrangian 

We construct the chiral Lagrangian using a spurion analysis of the quark level Lagrangian 
given in Eqs. (ITUi) and (ITTIl . We generally assume a power counting 

rUg ~ aA^ , (19) 

but work to the leading order necessary to parameterize the dominant artifacts from the 
anisotropy. At LO, the meson Lagrangian is given by^ 

4 = Ar (9^S9^St) - ^tr (m.BSt + S(m,B)t) 

£2 £2 

- ^—ti (a,WS^ + ^(a^W)^) - ytr (a^W^S^ + S(a,W«)'^) . (20) 
By taking functional derivatives with respect to the spurion fields in both the quark level 



and chiral level actions, one can show 37 



and similarly the new dimension-full chiral symmetry breaking parameters are defined as. 



W = Jim^ csw ^^ '""^'"^' , (22) 



f 



= lim ci^4-^^^^#^- (23) 

This anisotropic Lagrangian, Eq. (!20|) . has one more operator than in the isotropic limit 0,0], 
which is simply a reflection that there are now two distinct 0{a) operators in the Symanzik 
action. However, for a fixed anisotropy, ^ = as/ at, these two 0{a) operators are indistin- 
guishable and can be combined into one. Inserting the tree level values of the coefficients in 



the Symanzik action j29j, one finds 



W occsw = J^, (24) 

1 /o^ 

2 [at 



^ We remind the reader we are working in Euclidean spacetime. We are using the convention / ~ 132 MeV. 
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where the speed-of-hght parameter, u, is determined in the tuning of the anisotropic action. 
The first clear signal of the anisotropy begins at the next order, 0{ap'^) for the unimproved 
action and C(a^p^) for the improved action. We first discuss the unimproved case. 

These next set of operators are only present for the unimproved action. In the isotropic 
limit, it was shown there are five additional operators at this order 

C^,am = 2W^4tr tr (a.WSt + S(a,W)t) 

+ 2W^ii (a^SS^S^ [a,WE^ + S(a,W)^]) 
+ AW(iti (m^BS^ + S(m,B)1') tr (a.WST + S(a,W)1') 
+ AW-jii (m^BS^ - S(mgB)^) tr (a,WS^ - S(a,W)^) 
+ AWsti [niqWL^ as'WT) + S(mgB)^S(a,W)^) . (26) 

For the anisotropic action, there are an additional five operators similar to the above five 
with the replacement of low-energy constants (LECs) W^s ~^ W^l-s ^"^^ ^ simultaneous 
replacement of the condensate W —>■ W^. These five new operators, as with the 0{a) 
operators, are indistinguishable from those in Eq. fl26|) at a fixed anisotropy. There are two 
additional operators at this order however, which introduce new effects associated with the 
anisotropy, 

4,am = W^i^tr (5,S9.St) tr 

+ W^tr ((9^S(9^S^ u^ui [a.W^St + S(a,W«)1'] ) . (27) 

When the anisotropic vectors are set to their constant value, (u^)^ = (1,0), one sees that 
these operators lead to a modification of the pion (meson) dispersion relation. 



iEl+pl)[l + W^ + W^^j+ElW^^, (28) 
where 

W = 32{NfW4 + W5) , = 32{NfW^ + W^) , 

= lQ{NfW^ + Wl) , (29) 

and Nf is the number of fermion flavors. With the 0{a) improved anisotropic action, these 
effects all vanish and the leading lattice artifacts begin at 0{a'^). The chiral Lagrangian at 
this next order in the isotropic limit was determined in Ref. j3] for which there were three 
new operators. In the anisotropic theory, there are an additional six operators but just as 
with the 0{a) Lagrangian, the three new anisotropic operators can not be distinguished 
from those which survive the isotropic limit unless multiple values of the anisotropy are 
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used. The Lagrangian at this order is 



<t>,a 



AW, 



,Wtr (S + S^) 



ml 



1 2 



- m!j Lwtr (s - s^") 1 - 4iy| Lw« tr (S - s^) 

- 41^8 (a.W)^tr (SS + S^^S^) - AWl (a,W^)2tr (SE + T.^T.^) 

- AWl (a,W)(a,W«) [tr (S + S^) 1 ^ - 41^7^ (asW)(a,W^) [tr (S - T)) 

- AWl (a,W)(a,W«) tr (SS + S^E^) . 



(30) 



The last three operators in this Lagrangian vanish for an (9(a)-improved action as they 
are directly proportional to the product of the two 0{a) terms in the Symanzik Lagrangian, 
Eq. (fTOl) . The other six operators in this Lagrangian receive contributions both from products 
of the 0{a) Symanzik operators as well as from terms in the 0{a?) Symanzik Lagrangian.^ 
Therefore they are still present for an (9(a)-improved action but the numerical values of 
their LECs, W^_^ and W^|_8, will be different in the improved case. 

At C(aV), there are nine new operators, three of which survive the isotropic limit, and 
three of which explicitly break the accidental 0(4) symmetry down to 0(3), 



Qi(a,W)2tr (a^Sa^St) + Q^ia^^mfti {d^Ld^L^) tr (S + S^) 

+ g«(a, W«)2tr (9^S9^St) tr (S + S^) + g|(a,W«)'tr (S^SS^S^ [S + S^] ) 
+ Qiias'm^fu^.ui tr {d^T.d,T)) + Qi{a,m^fu^^ul tr {d^T.d,T)) tr (S + S^) 



+ Q|(a,W«)\jM« tr (9^S9,St [S + T)]) 



(31) 



The first operator in each set of three with coefficients Qi, Q\ and Q\ are modifications of 
the LO kinetic operator, while the remaining operators additionally break chiral symmetry. 
The first operator in Eq. fl3Tl) is an example of the operators mentioned before which do not 
break any of the (approximate) lattice symmetries, and are therefore amount to polynomial 
renormalizations of the continuum LECs. In this example, we see with the replacement 



f + 8gi(a,W)^ 



(32) 



the entire effects from this 0{a?p'^) operator are renormalized away to all orders in the EFT. 
In the anisotropic theory, this also works for the operator with LEG Qi, as this operator 
does not explicitly break 0(4) symmetry. This provides a further example of effects which 
arise because of the anisotropy but contribute to the low-energy dynamics in an isotropic 
fashion. The operators with explicit anisotropic vectors will modify the dispersion relation 
as in Eq. (!28|) but at 0{a?). For the C(a)-improved action, these operators will provide the 
dominant modification to the pseudo-Goldstone dispersion relations. 

The last set of operators we wish to address for the meson Lagrangian are those which 
explicitly break the accidental 0(3) symmetry down to the cubic group. There are two 



^ For conventional reasons 1^] we have normalized these operators to the square of the condensates appearing 
at 0{a), but one should not confuse this to mean that these operators vanish for the C'(a)-improved action. 
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operators in the meson chiral Lagrangian but they stem from only one quark level operator 
at 0{a^) and therefore all the LO 0(3) breaking effects for all hadrons can be removed with 
the inclusion and tuning of one new operator in the action, ^04^'' from Eq. f|T7j) . The 0{3) 
breaking operators in the meson chiral Lagrangian are 

^lo(3) = CKasW^)%AJlA tr (9,,Sa.St) tr (5,S9.St) 

+ Clia,Wi)%JlSl5l tr (a,,Ea.Eta,Ea.Et) , (33) 



1. Aoki Regime 

Aoki first pointed out the possibility that at finite lattice spacing, lattice actions can 
undergo spontaneous symmetry breaking of flavor and parity in certain regions of phase 
space 131 . In Ref. 0], Sharpe and Singleton addressed this possibility within an EFT 



framework by extending the meson chiral Lagrangian to include lattice spacing contributions. 
We begin by summarizing the discussion of Sharpe and Singleton which will allow us to 
highlight the new effects that arise from the anisotropy. For clarity of discussion, we consider 
the unimproved two-flavor theory in the isospin limit. Following the notation of Ref. 0], the 
non-kinetic part of the chiral potential can be written 

Vx = -|tr(E + Et) + |[tr(E + Et)]^ (34) 

where ci and C2 are functions of the quark mass nig and the lattice spacing a, 

rriq 
A 

-j^ + mga + a^A^\+ ... . (35) 

The ". . . " denote higher order terms in the quark mass and lattice spacing and we are ignor- 
ing dimensionless numbers of 0(1). Assuming a power counting niq ~ aA^, the contributions 
to the vacuum from the C2 term are suppressed compared to the Ci term, and the vacuum is 
in the continuum phase with the chiral condensate aligned with unity. The Aoki-phase can 
occur when there is a fine-tuning of the quark mass and lattice spacing contributions to Ci 
such that overall size of ci becomes comparable to C2. Parameterizing the S-field as 

S = A + 2r • B , (36) 

with the constraint + = 1, the potential, Eq. flMl) becomes 

= -ciA + C2A^ . (37) 

If the minimum of the potential occurs for — 1 < < 1, then the vacuum 

So = ( S ) = Ao + zr ■ Bo , (38) 

develops a non-zero value of Bq, spontaneously breaking both parity and the remnant vector- 
chiral symmetry, SU{2)v — * U{1) [2], giving rise to one massive pseudo-Goldstone pion and 
two massless Goldstone pions. If the minimum of the potential occurs for .A^j^j ^ — 1 or 
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Amin > 1 then the vacuum hes along (or opposite) the identitiy, |y4o| = 1, |Bo| = with the 
same symmetry breaking pattern as QCD. 

For unimproved Wilson fermions in the isotropic limit, the leading contributions to ci are 

ci = f (mqM + a,w) . (39) 

In the anisotropic theory, there is an additional contribution to the LO potential, such that 

ci cf = ci + fasW^ 

= f (mqM + a,W + a,W^) . (40) 

If the two terms which contribute to W oc 2{qatiFtiq + qaijFijq) are of opposite sign, then 
oc {qauFuq) may be the dominant lattice spacing contribution to Ci. Therefore, the 
anisotropic theory may be in the Aoki-regime even when the isotropic limit of the theory 
(with the same a^) is not, and vice versa. This same discussion holds for (9(a)-improved 
actions as well but with a different power counting. 

This analysis carries important consequences for anisotropic actions with domain-wall 
and overlap fermions as well. These actions are generally tuned to lie between the first two- 
fingers of the Aoki-regime in the mg — plane |38| . The optimal value of the bare fermion 
mass, which provides the least amount of residual chiral symmetry breaking, may be shifted 
in the anisotropic theory relative to the isotropic value, and furthermore the allowed values 
of the coupling for which there is a QCD-phase may shift as well. 



B. Heavy Baryon Lagrangian 



We now construct the operators in the baryon Lagrangian which encode the leading lattice 

and explicitly 



artifacts from the anisotropy. We use the heavy baryon formalism [3^ 
construct the two-flavor theory in the isospin limit including nucleons, delta-resonances and 
pions. The extension of this to include the octet and decuplet baryons is a straight forward 
exercise. This construction builds upon previous work in which the heavy baryon Lagrangian 
has been extended in the isotropic limit to include the 0{a) ^ and C(a^) ^ lattice artifacts 
for various baryon observables.^ The Lagrangian is constructed as a perturbative expansion 
about the static limit of the baryon, treating it as a heavy, static source as with heavy 



quark effective theory [i^, 3]. In building the heavy baryon chiral Lagrangian it is useful 



to introduce a new chiral field, 



a 



exp 



(41) 



In those works @] , the heavy baryon Lagrangian was extended exphcitly for Wilson fermions and in the 
latter also for a mixed action scheme with chirally symmetric valence fermions. In this work, we explicitly 
address the anisotropic effects for these two scenarios. One can also include the anisotropic effects for 
twisted mass 4l| and staggered fermions 4^ |43| following the heavy baryon construction in Refs. 44 1 



and 



45 



46[ respectively. However, we are unaware of specific plans to generate anisotropic twisted mass 



or staggered fermions so we do not pursue this here. 
® This field is generally denoted as ^ but to avoid confusion with the anisotropy parameter, we use a. 
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which transforms under chiral rotations as 



LaV^ = VaR^ , 



(42) 



with V defined by Eq. fH21) . One then dresses the baryon fields with cr such that under chiral 
transformations, the nucleon and delta fields transform as 



ViitVjjiVkk'Ti'jik' , 



(43) 



where is a two-component fiavor field and T is a fiavor-symmetric rank-three tensor, 
normalized such that 
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A 



++ 



n, 



Ti22 = ^A' 



222 



A" 



(44) 



Both the nucleon and delta are treated as heavy matter fields and are constrained by the 
relations 



N 



T 



2 
2 



(45) 



where is the four velocity of the baryon which can be chosen in its rest frame to be 
Vfj^ = (1, 0). This has the effect of projecting onto the particle component of the field in the 
rest frame of the baryon. The spin-3/2 fields can be described with a Rarita-Schwinger field, 
which in the heavy baryon formalism gives rise to the constraints 



v-T = 



S-T = 



(46) 



where 5*^ is a covariant spin vector [39|, |40| . The simultaneous inclusion of the nucleon and 
delta fields introduces a new parameter into the Lagrangian, the delta-nucleon mass splitting 
in the chiral limit, 



rriT — mN 



(47) 



which is generally counted as A ~ in the chiral power counting [39|, |40|, |49| . Because 



this mass parameter is a chiral singlet, it leads to a modification of all the LECs in the 
heavy baryon Lagrangian. These particular effects can be systematically accounted for be 
treating all LECs as polynomials in A [s^, 51, 5^, 53|. In the isospin limit, the heavy baryon 
Lagrangian, including the leading quark mass terms, the leading lattice spacing terms and 
the leading baryon-pion couplings is given by 

= NvDN + atasNiv ■ u^u^ ■ D)N - 2aMNmi{M+) 

+ f^[v D + /\]T^ + atCXET^iv ■ u^u^ ■ D)T^ - 2aM%T^ti{M+) 
+ 2gANS ■ AN - 2g^^f,S ■ AT, + (^ajv \TfA^^,,Nk + h.c] . 
- 2awNmi{W+) - 2ai,Nmi{Wi) - 2a^^7;T^tr(W+) - 2a^f^r^tr(>Vi) (48) 
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(a) (b) (c) (d) 

FIG. 1: Diagrams contributing to the nucleon mass at LO ((a) and (h)) and NLO ((c) and (d)). 
Figure (a) is an insertion of the leading quark mass term proportional to gm- Figure (h) is an 
insertion of the lattice spacing terms proportional to ay/ and a^. The loop graphs arise from the 
pion-nucleon and pion-nucleon- delta couplings. These loop graphs generically scale as but also 
depend upon A, and away from the continuum limit upon the lattice spacing as well. All vertices 
in these graphs are from the Lagrangian in Eq. 

In this Lagrangian, the chiral symmetry breaking spurions are given by 

M+ = ^ (o-mjcr + a^mgcx^) , 
n;+ = y ((tWV + a^Wa^) , 

= y ((t(W«)V + (j^WV^) , (49) 

the chiral covariant derivative is 

(D^T),,, = d^T,,u + V^,%fc + Vj;.,T,,v, + V,^,,T,,fc, (50) 
and the vector and axial fields are respectively given by 

A,= '- {ad^a^ - a^d^a) . (51) 

Relative to the Wilson extension of the heavy baryon Lagrangian in the isotropic limit [1, 0] 
there is one additional mass operator for the nucleon and delta fields. There are additionally 
extra derivative operators which give rise to the leading modification of the dispersion rela- 
tion for the baryon fields. As in the meson chiral Lagrangian, Eq. fl2Ul) . at fixed anisotropy, 
these two additional mass operators (with coefficients cr^ and a^) can not be distinguished 
from their counterparts which survive the isotropic limit (with coefficients aw and aw)- At 
this order in the lattice spacing, 0{a), the baryon-pion couplings, qa, Qan and qaa, are not 
modified This Lagrangian, Eq. (148|) gives rise to the LO and NLO mass corrections to 
the baryon masses. For example, in Fig. [1] we display the graphs contributing to the nu- 
cleon mass at 0{mg), Fig. [I](a), 0{a), Fig. [I](b) and 0{mg^'^), Fig. W^c) and Fig. W^d). The 
expressions for the quark mass dependence of the nucleon and delta masses can be found 
for example in Ref. [52|, and the lattice spacing dependent corrections to these masses for 
Wilson fermions in the isotropic limit can be found in Refs. jB], 0] . 

At 0{amq), there are two types of mass corrections, those from loop graphs from the 
lattice spacing dependent operators in Eq. fHSj) and tree level terms from operators in the 
0{amq) Lagrangian. These corrections will be identical in form to those present in the 
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isotropic limit, but with different numerical values due to the new anisotropic operators 
which do not explicitly break the hypercubic symmetry. The loop graphs for the nucleon 
and delta mass corrections are depicted in Ref. [6i] in Figs. 1 and 2 respectively. In the 
isospin limit, there are two tree level operators for the nucleon and deltas each, 

J^ntI = - nMwNmi{M+)ii{W+) - nij^Nmi{M+)ii{yvi) 

+ tMH^f^T^tr(A^+)tr(W+) + 4^^f^T^tr(A^+)tr(Wi) , (52) 

which also function as counterterms for divergences from the above mentioned loop graphs. 
As with the mesons, we must consider the 0{a?) effects to see the first explicit breaking of 
the anisotropy. Even in the isotropic limit, the lattice spacing corrections to the baryons at 
this order proved to be quite interesting, with the presence of the first operators which arise 
from the 0(4) breaking down to the hypercubic group at this order [6|. There are several 
generic 0{a?) operators which do not explicitly break the hypercubic group, 

i^ir^ = - nwwNmi{W+)ti{W+) - n\^^N mi{Wi)ti{Wi) 
+ tH/w'T^T^tr(>V+)tr(W+) + t\^^%T 

- n\,^Nmi{W+)ti{Wi) + 4^^f^T,,tr(W+)tr(Wi) , (53) 

where the last two operators vanish for the 0{a) improved action and the first two sur- 
vive the isotropic limit. The 0(4) breaking operators appear at this order in the heavy 
baryon Lagrangian because in the heavy baryon theory, there is an additional vector, the 
four-velocity of the baryon one can use to construct operators. For the mesons, the only 
vector is dn and so the 0{A) breaking operators in the meson Lagrangian don't appear until 
0{a'^p^) [4'|. These heavy baryon 0(4) breaking operators are j6| 

^ntI = Nv^v^Vf,Vf,N + a% fpV^v^v^v^jTp . (54) 
Similar to these, the 0(3) breaking operators appear at this order, 

^nt\ = «'^3 NSlJlJlZl N + a% %5l5l5l5l T, . (55) 

However, even though these two sets of operators explicitly break the 0(4) and 0(3) sym- 
metries respectively, they do not lead to modifications of the dispersion relations until higher 
orders and thus function as mass corrections. There is a subtlety of the heavy baryon 
Lagrangian related to reparameterization invariance (RPI) jH^], which constrains coefficients 
of certain operators in the heavy baryon Lagrangian. For example, in the continuum limit, 
the LO kinetic operator of Eq. (HHI) is related to a higher dimensional operator in such a 
way as to provide the correct dispersion relation,^ 

n2 

C = Nv-DN — ^Nv-DN + N—^N , (56) 
such that the energy of the non-relativistic nucleon is given by 



In this equation, D^^ — D"^ - {v ■ DY . 



14 



For the anisotropic theory, this dispersion relation is then modified by an additional operator. 
While in the continuum limit, RPI constrains the coefficient in front of the 1/Mn operator, 
the anisotropic action gives rise to modifications of this relation. 

C = NvDN + N—^N — > 
2Mn 

_ n2 

= Nv DN + N—^N + atasNv ■ u^u^ ■ DN + as(JKEN{u^ ■ DfN . (58) 
The resulting nucleon dispersion relation from this action yields 

En = Mm • ^ 



2Mn 

En{1 + atas) = M^(l + atas) - 2Nfas{aw + afj.) + + a.aKE {En - M^f • (59) 

For the 0(a)-improved action, the form of this dispersion relation stays the same with 
{at, as} — > {al^al}. Comparing this dispersion relation to that of the pion, Eq. (128|1 . it is 
clear that even if the pion dispersion relation were tuned to be continuum like, the nucleon 
dispersion relation would still contain lattice spacing artifacts. This is a simple consequence 
of the fact that the LECs in the heavy baryon Lagrangian (both the physical and unphysical) 
have no relation to those in the meson chiral Lagrangian. 



IV. DISCUSSION 

We have developed the effective theory describing pions, nucleons and deltas at low 
energy, including the finite lattice spacings effects from an anisotropic Wilson lattice action. 
In particular, we focussed on the correction to the dispersion relations caused by the breaking 
of hypercubic symmetry. The theory considered here is suited for the anisotropic Wilson 
lattices currently being generated. Extensions to a mixed action case [s^, 5^, for example 



with domain-wall or overlap valence fermions and the Wilson sea fermions, both anisotropic, 
can be made with ease as in Refs. HHHIej. We have also highlighted some subtlety 



involving the Aoki-phase for anisotropic actions, which may be of importance in tuning the 
bare fermion masses in domain-wall and overlap fermions; for a fixed bare fermion mass and 
spatial lattice spacing, if the isotropic action is in the QCD-phase, the anisotropic theory may 
still be in the Aoki-phase. This work will aid in the continuum extrapolation of correlation 
functions computed with anisotropic lattices. 
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